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ABSTRACT 

We have acquired simultaneous high-precision space photometry and radial 
velocities of the bright hybrid (3 Cep/SPB pulsator 7 Peg. Frequency analyses 
reveal the presence of six g modes of high radial order together with eight low- 
order 13 Cep oscillations in both data sets. Mode identification shows that all 
pulsations have spherical degrees I = — 2. An 8.5 Mq model reproduces the 
observed pulsation frequencies; all theoretically predicted modes are detected. 
We suggest, contrary to previous authors, that 7 Peg is a single star; the claimed 
orbital variations are due to g-mode pulsation. 7 Peg is the first hybrid pulsator 
for which a sufficiently large number of high-order g modes and low order p and 
mixed modes have been detected and identified to be usable for in-depth seismic 
modeling. 

Subject headings: stars: oscillations - stars: variables: other - binaries: spectroscopic 
- stars: early-type - stars: individual (7 Peg) 
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Introduction 



The bright {V = 2.8) B2 IV star 7 Peg was recog nized as a pulsating variable of the 



j3 Cep class more than 50 years ago (McNamara 



periodic until recently. 



Chapellier et al 



19531 ). and it was believed to be singly 



(I2OO6I ) studied the star spectroscopically and 
demonstrated its multiperiodicity. These authors also examined the claim that 7 Peg is 
a spectroscopic binary a nd deduced an eccentric (e = 0.6 2) 370.5-day orbit. This orbital 



solution was disputed by 



Butkovskaya fc Plachindal (120071 ) who suggested that the orbital 



eccentricity was spurious and caused by outbursts similar to those of Be st ars, and favored 



an orb ital period near 6.8 d, in accordance with the original suggestion by 



f!1979h . 



Harmanec et al 



The multiperiodic oscillations of 7 Peg are highly interesting beca use they are caus ed 



by tw o different sets of pulsation modes: two frequencies detected by 



Chapellier et al. 



(120061 ) correspond to low-order pressure (p) and gravity (g) modes typical for (3 Cep 
stars, but the other two are high-order g modes as excited in the Slowly Pulsating B 
(SPB) stars. Indeed, 7 Peg is lo c ated i n the overlap region of both types of variables 



in the HR diagram (see iHandlerl (120091 )). The frequencies of the two sets of modes 
are sensitive to the physical conditions in different parts of the stellar interior. Hybrid 
oscillators therefore offer the possibility to obtain a more complete picture of the physics 
inside a star using asteroseis r nology where pulsations act as seismic waves (see, e.g., 
Dziembowski fc Pamyatnykhl (120081 ) for case studies). 



However, the possible binarity of 7 Peg imposes difficulties. The detection of low- 
frequency oscillations can be compromised by an inaccurate orbital solution. Photometry 
does not suffer from this problem (provided the light-time effect is negligible). Consequently, 
Handler! (120091 ) carried out a multicolor time-series photometric study of 7 Peg, detected 
four SPB-type pulsation modes, and confirmed the two modes of /3 Cep type. One of the 
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latter was identified as radial; it would be either the fundamental mode or the first overtone, 
immediately constraining the mean stellar density. The enormous asteroseismic promise of 
7 Peg motivated us to perform a high-precision photometric and radial velocity study of its 
pulsations, in the hope to detect a sufficient number of p and g modes to sound its interior 
structure. 



2. Observations 



7 Peg was observed with the MOST satellite (jWalker et al.ll2003l ) from Sep 16 to Oct 



16, 2008. Due to the brightness of the target. 



reduction method developed by 



Reegen et al. 



'abry imaging mode was used, and the data 



(I2OO6I ) to minimize the effects of stray light 



was employed. There were over 55000 data points at a cadence of 30 s in the original data 
set. These were summed into 4018 0.005-d bins having an rms scatter of 1 mmag per point 
and an effective Nyquist frequency of ~92 cd~^ . 

A simultaneous ground-based spectroscopic multisite campaign was organized for 
pulsational mode identification. The majority of the data originate d from the Automati c 



Eaton &: Williamson 



Spect roscopic Telescope of Tennessee State University (TSU-AST, 
20071 ). and consisted of 1660 spectra (4900-7100 A) taken over a span of ten weeks (Sep 3 
- Nov 15, 2008) with an effective Nyquist frequency of ~68cd^^ . We have reduced these 
spectra and extracted radi al velocities from 34 lines with the techniques described by 
Eaton fc Williamsoru (120071 ). The external error of these radial velocities is about 0.2 kms^^ . 



This high precision results from 7 Peg having a very sharp lined spectrum for a hot star 



biting et al 



(I2OO6I ) listed V sin z = kms ^ ! Concer ning the star's metallicity, 



( 2006 ) derived Z = 0.0091 ± 0.0021 from optical and 



Morel et al. 



Daszyhska-Daszkiewicz fc Niemczura 



(120051 ) determined [m/H] = —0.04 ih 0.08 from ultraviolet spectra. 



- 7- 



Additional spectroscopy was carried out at three observatories, and ground-based 
time-resolved multicolor photometry was acquired at three more sites. However, the present 
paper only reports the initial results from the MOST space photometry and TSU-AST 
radial velocities. 



3. Frequency analysis 



The h eliocentrically correcte d data were searched for periodicities using the program 



Period04 (ILenz fc Breger 



20051 ). Amplitude spectra were computed, compared with the 
spectral window functions, and the frequencies of the intrinsic and statistically significant 
peaks in the Fourier spectra were determined. Multifrequency fits with all detected signals 
were calculated step by step, the corresponding frequencies, amplitudes and phases were 
optimized and subtracted from the data before computing residual amplitude spectra, 
which were then examined in the same way. 

This analysis was performed for the MOST and radial velocity data independently. We 
conservatively only accepted signals that exceeded an amplitude signal-to-noise ratio of five 
in at least one of the data sets and that were prominent in the other. Some steps of this 
procedure are shown in Fig. [H Fourteen independent signals were detected; the agreement 
between the photometric and radial velocity measurements is remarkable. The residual 
amplitude spectrum after this solution was featureless for the radial velocities. Some peaks 
in the residual MOST data remained, but to err on the side of caution we did not consider 
them to be intrinsic to 7 Peg. 

The frequencies of all these signals are consistent within the errors between the two 
data sets. We have d etermined weighted mean values of t hose frequencies by computing the 



formal uncertainties (iMontgomery &: O'Donoghue 



I999I ) in the individual data sets and 
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then applying their inverse squared as the weight. With the resulting frequencies fixed, we 
have re-calculated the amplitudes, phases and amplitude signal-to-noise ratios and list the 
results in Table [TJ 

All but the two closest frequencies are resolved within our data set. Their difference 
/i2 — /s = 0.0112 cd^Ms 82% of the time resolution of the radial velocity data. Since 
a significant peak remains after prewhitening /s we accept fu, but caution that the 
parameters of these two variations in Table [1] may have systematic errors. 



4. Discussion 



We start by examining previous claims that 7 Peg is a spectroscopic binary. Our radia l 



veloci ties were acquired in a period of time where the orbital solution by 



Chapellier et al. 



(120061 ) predicts a change in radial velocity of about 20 kms ^ . Our ri i ghtly mean radial 



velocities are constant within ±0.3 kms ^ 



Butkovskaya fc Plachindal (120071 ) claimed an 



orbit with a 6.816 d period and about 0.8 kms^^ radial velocity amplitude. Again, this is 
inconsistent with our data. Instead, the 1 cd~^ alias of this "orbital" period lies well within 
the domain of the g-mode frequencies (1 — 1/6.816 d=0. 8533 cd~^ ); the strongest g-mode 
pulsations have radial velocity amplitudes of 0.7kms~^ . We conclude that 7 Peg is not a 
6.8-d spectroscopic binary. The cause of the spora dic radial velocity changes rema ins to be 



understood. The possibility of Be-star outbursts (IButkovskaya fc Plachinda 



20071 ) seems 



remot e given that 7 Peg likely rotates intrinsically slowly, unlike the Be stars (e.g.. 



Porter 



(U996|)). 



According to the position of 7 Peg in the HR diagram determined by 



Handlei] (120091 ). 



all signals with frequencies below led ^ are due to high -order g r aodes ; the remaining 



variations are caused by low-order p and mixed modes. 



Handleil (|2009|) identified /i as 
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a radial mode and argued that /s is a dipole mode; these are the only two modes with 
reasonably secure identifications. From comparison with pulsation models this author 
concluded that models in just two small domains of parameter space in the HR diagram 
are consistent with these mode identifications and with the frequency values. In case /i 
is the radial fundamental mode, models around 8.5 M© match the observations. If /i 
corresponded to the first radial overtone, models with masses a round 9.6 Mq reproduce /i 



and /s. All models in the present paper and by 



Handler! (120091 ) were computed with the 



the W arsaw-New Jersey stellar evolution and pul sation code (e.g . , see 



Asplund et al. 



Pamyatnykh et al. 



( 120041 ) element mixture. 



(119981 ) for a description), using OP opacities, the 
and providing linear nonadiabatic model frequencies. No rotation or convective core 
overshooting was included in the models for simplicity of this preliminary model fitting. 
Figure [2] compares the additional i = — 2 P Cep-type pulsation frequencies predicted by 
the two models with our new observations. 

All newly detected pulsation frequencies are explained by the 8.5 Mq model (upper 
panel of Fig. [2]). At first sight, there are two small inconsistencies: the close doublet 
/5//12 only has one theoretical counterpart, and /13 is not well matched by the theoretical 
value. The doublet can be explained by possible rotational splitting at Vrot ~ 3 kms~^ , 
consistent with the very low v sin i of the star. The mismatch for /13 may be more 
of a p roblem, but it would not be new: in their an alysis of the pulsation spectrum of 



u Eri, iPamyatnykh. Handler, fc Dziembowskil (|2004| ) noted the same problem for the 
highest-frequency p mode {£ = 1, P2) observed. Interestingly, in the present model /13 
would also correspond to the £ = 1, p2 mode and even the size of the frequency mismatch is 
similar. It can be suspected that this mismatch originates from inadequate physics in the 
models that can be improved through asteroseismology. 



The 9.6 Mq model (lower panel of Fig. [2]) results in poorer agreement between the 
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observed and theoretical frequencies, although the number of theoretically predicted modes 
in the domain of excited frequencies is larger. In particular, only one theoretical i = 1 mode 
is available to match fiQ and /14. One would need to invoke rather fast rotation to explain 
both modes, but in this case 7 Peg would be viewed close to pole-on, which would cause 
heavy geometrical cancellation of just these pulsations. We conclude that, in all likelihood, 
7 Peg is a ~ 8.5 Mq star oscillating with a dominant radial fundamental mode. In this 
case, all the theoretically predicted i = — 2 j3 Cep-type pulsation modes are observed. 

The amplitudes and phases of the oscillations can be used to type the pulsation modes. 
We have computed theoretical radial velocity to light amplitude ratios and phase differences 
with for the M O^'T bandpass. The 8.5 M0 model was used together with static atmospheres 



( iKuruczl l2004l ) with a metallicity parameter [m/H] = 0.0 and a microturbulence velocity 
= 2 km/s. The comparison of the theoretical and observed amplitude ratios and phase 
shifts is shown in Fig. [3l modes with < i < 4 were considered. 

The measured amplitude ratios and phase shifts are consistent with the interpretation 
that all oscillation frequencies are caused by modes with £ < 2; theoretical results for modes 
with higher i are off scale in Fig. [31 This corroborates the identification obtained in Fig. [2j 

Together with our ground-based multicolor photometry, our radial velocities can 
be used to constrain nonadiabatic pulsation theory, as empirical determinations of the 
/ parameter used to describe the bolometric flux amplitude depending on the surface 
displacement can be m ade. In turn, this can be used to choose the most suitable o pacities 



for model calculations (jPaszyhska-Daszkiewicz. Pamyatnykh. fc Dziembowski 



20051) 



Figure H] shows a comparison of theoretically predicted g-mode frequencies of our 
8.5 M0 model with the observations. The frequency separations of the g modes are well 
explained by this model: five of the six observed modes may form a sequence of consecutive 
radial overtones of £ = 1 modes. The sixth mode would then he i = 2. The i = 1 modes 
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are stable; a final seismic model must explain their excitation. 

We conclude that 7 Peg presents sufficient information to carry out detailed seismic 
modehng of a hybrid pulsator for the first time. The eight P Cephei pulsation frequencies 
restrict the possible models and their parameters considerably. The observed amplitude 
ratios and phase shifts between the radial velocities and photometric data provide clues 
towards the opacities to be used. Theoretical models constrained by such a large set of 
observables must also reproduce the high-order g modes and the excited frequency domains. 
The present study is a demonstration of the value of combining space photometry with 
ground-based spectroscopy of bright stars. 

This work has been supported by the Austrian Fonds zur Forderung der wis- 
senschaftlichen Forschung under grant P20526-N16. AFJM is grateful for financial support 
from NSERC (Canada) and FQRNT (Quebec). We thank Frank Fekel of Tennessee State 
University for providing the list of photospheric lines we used for extracting velocities 
from the TSU spectra. We are grateful to Ramotholo Sefako for organizing the SAAO 
observations. 
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Fig. 1. — Amplitude spectra of our 7 Peg data with consecutive prewhitening. The pho- 
tometric data are compared with the radial velocity measurements (reverted graphs). The 
frequency regions 1.5 — 5.5 and > 10 cd^^ are not shown as they c ontain no intrinsic signals. 



Handler] t00% . The ap 



The frequencies /i — /g have been assigned for consistency with 
parent 1 cd~^ sidelobe of /i in the MOST data is an alias peak originating from the orbital 
frequency. 
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Fig. 2. — Comparison of theoretically predicted pulsation frequencies of two models fitting /i 
and /s. The model frequencies are plotted at the bottom of the panels, with amplitudes of 0.7 
(radial modes), 0.5 (dipole modes), and 0.33 (quadrupole modes). The observed frequencies 
protrude from the top of the panels. 
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Fig. 3. — Observed (full circles with error bars) and theoretically predicted radial-velocity 
to photometric amplitude ratios and phase shifts for the fourteen signals detected. The 
upper panel shows the /3 Cep pulsations and the lower panel the SPB oscillation modes. The 
star symbols show the theoretical locations of radial modes in this diagram, the diamonds 
represent dipole modes and the triangles stand for quadrupole modes. Only model modes 
with frequencies in the observed domains are plotted. 
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Fig. 4. — Comparison of tiicorctically predicted g-mode pulsation frequencies for a 8.5 M© 
model of 7 Peg with the observed frequencies. As ordinate the stabihty parameter 77 is used; 
if it is larger than zero, the corresponding mode is excited in the model. Diamonds represent 
dipole modes and triangles stand for quadrupole modes. The vertical lines are drawn at and 
identified with the observed frequencies. 
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Table 1. Multifrequency solution for our photometric and radial velocity data. 



ID Frequency Amplitude Amplitude {S/N)phot {S/N)rv 4>Vrad - <Pmost 



(cd ^ ) (mmag) (kms ^ ) (rad) 



/l 


6.58974 ± 0.00002 


6.59 


3.359 


300.9 


327.9 


4.318 


± 


0.004 


/2 


0.63551 ± 0.00010 


1.70 


0.500 


21.3 


17.5 


2.189 


± 


0.019 


h 


0.68241 ± 0.00007 


1.99 


0.736 


24.8 


26.0 


2.280 


± 


0.015 


h 


0.73940 ± 0.00010 


1.23 


0.522 


15.4 


18.6 


2.481 


± 


0.023 


h 


6.01616 ± 0.00014 


1.14 


0.358 


45.0 


34.4 


4.55 


± 


0.03 


h 


0.88550 ± 0.00007 


0.90 


0.723 


11.3 


25.9 


3.18 


± 


0.03 




6.9776 ± 0.0005 


0.33 


0.095 


16.3 


9.5 


4.40 


± 


0.10 


h 


0.91442 ± 0.00011 


0.51 


0.464 


6.4 


16.6 


3.79 


± 


0.05 


h 


6.5150 ± 0.0008 


0.21 


0.063 


9.4 


6.3 


4.78 


± 


0.15 


fio 


8.1861 ± 0.0008 


0.18 


0.064 


9.2 


6.7 


4.16 


± 


0.17 


hi 


0.8352 ± 0.0003 


0.27 


0.180 


3.4 


6.3 


2.61 


± 


0.09 


/l2 


6.0273 ± 0.0005: 


0.33 


0.112 


12.4 


10.3 


4.81 


± 


0.10 


/l3 


9.1092 ± 0.0012 


0.12 


0.041 


5.9 


4.5 


3.97 


± 


0.26 


/l4 


8.552 ± 0.002 


0.10 


0.027 


5.0 


3.2 


4.06 


± 


0.34 



Note. — The formal errors on the photometric amplitudes are ± 0.02 mmag; the 1-a errors on 
the radial velocity amphtudes are ± 0.007 kms~^ . /12 has been marked with a colon because it 
is not fully resolved from f^. 



